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1.  Introduction 

The  millimeter- wave  spectrum  holds  promise  for  a  wide  range  of  existing  and  proposed 
military  and  civilian  electronics  systems.  Satellite  and  terrestrial  communications,  remote 
sensing,  and  weather  and  air-traffic  control  radars  all  benefit  from  the  small  antennas,  wide 
bandwidths,  and  high  imaging  resolution  offered  by  millimeter-wave  systems.  Realization 
of  the  full  potential  of  many  such  systems,  however,  awaits  the  development  of  practical 
and  reliable  high-power  solid-state  sources  and  amplifiers.  At  present,  the  generation  and 
amplification  of  high-power  mii  :eter-wave  signals  is  dominated  by  vacuum-tube  devices 
such  as  Gyrotrons  and  Travelling  Wave  Tubes  (TWTs)  (1).  This  is  undesirable  in  many 
applications  due  to  the  large  size  and  weight  of  the  systems,  and  the  high-voltage  power 
supplies  required  by  the  devices. 


Figure  1  —  Output  power  from  a  variety  millimeter-wave  sources.  This  figure  shows 
clearly  that  large  numbers  of  semiconductor  devices  arc  required  to  produce  powers 
comparable  with  tubes. 

In  order  to  compete  with  millimeter-wave  TWTs,  solid-state  transmitters  must  use  large 
numbers  of  devices  (figure  1).  Several  different  approaches  to  the  power-combining  of  multi¬ 
ple  millimeter-wave  devices  have  been  reported  (2).  Traditional  combining  techniques  using 
hybrid  3  dB  couplers,  or  large  numbers  of  devices  in  a  resonant  cavity  have  fundamental 
limits  regarding  efficiency  and  the  number  of  devices  that  can  be  combined.  These  limita- 
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tions  arise  from  unavoidable  circuit  losses  and/or  the  requirement  of  a  small  cavity  in  order 
to  avoid  multi-mode  problems.  Therefore  these  techniques  are  not  suitable  for  combining 
large  numbers  of  devices  at  millimeter-wave  frequencies. 

In  recent  years  a  planar  quasi-optical  approach  has  been  suggested  for  combining  the  output 
powers  of  millimeter-wave  sources  [3j.  In  this  approach  a  two-dimensional  array  of  semi¬ 
conductor  devices  is  used,  as  in  figure  2.  Each  device  is  integrated  with  a  printed  antenna, 
coupling  it  to  the  radiation  field.  In  this  way  the  power-combining  takes  place  in  free-space, 
and  thus  high  combining  efficiencies  (close  to  100%)  are  possible.  The  transverse  dimension 
of  quasi-optical  systems  can  be  quite  large,  which  accommodates  many  devices.  Devices 
need  not  be  tightly  packed  together,  so  removal  of  waste  heat  is  facilitated.  Furthermore, 
where  traditional  combining  techniques  require  many  non-reciprocal  elements  to  prevent  de¬ 
vice  interaction,  the  quasi-optical  arrays  actually  depend  on  the  controlled  interaction  of  the 
devices  for  proper  operation.  These  arrays  are  expected  to  have  application  in  a  wide  range 
of  frequencies,  up  to  the  TeraHertz  range  using  suitable  devices  [4]. 


and  printed  antennas  Ground  plane 

and  Heatsink 


Figure  2  —  Planar,  solid-state  millimeter-wave  array.  Large  output  powers  may  be 
obtained  at  millimeter-wave  frequencies  by  quasi-optically  combining  the  outputs  of  a 
large  number  of  semiconductor  devices. 

Prior  to  the  work  performed  in  this  grant  there  was  one  main  quasi-optical  array  demon¬ 
stration,  originating  with  Rutledge  at  the  California  Institute  of  Technology,  consisting  of  a 
distributed  oscillator  system  in  which  devices  are  mounted  in  a  periodic  grid  structure  [5-8] 
and  placed  in  an  open  quasi-optical  cavity.  It  has  been  pointed  out  [9]  that  this  approach 
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is  analagous  to  a  laser,  in  which  the  distributed  oscillator  system  acts  as  the  active  gain 
medium,  with  the  external  reflectors  providing  the  feedback.  A  100-element  MESFET  array 
has  been  demonstrated  successfully  using  this  approach  [10]. 

A  second  approach,  originating  from  the  work  reported  here,  involves  arrays  of  weakly- 
coupled  individual  oscillator  elements  [11].  This  system  forms  a  classical  antenna  array 
in  which  each  radiating  element  is  itself  a  free-running  oscillator.  The  array  elements  are 
synchronized  through  weak  mutual  coupling  mechanisms.  The  strength  of  the  coupling 
between  elements  is  limited  to  ensure  that  each  element  operates  close  to  its  free-running 
state,  hence  the  operating  frequency  is  solely  determined  by  the  design  of  the  individual 
oscillator  elements. 

Coupled-oscillator  arrays  may  have  advantages  over  the  distributed  oscillator  grids.  For 
one,  the  array  design  is  simplified  by  the  use  of  weak  coupling,  since  each  element  of  the 
array  is  designed  as  a  free-running  oscillator  and  is  only  slightly  perturbed  from  this  state 
during  operation.  The  technique  is  thus  modular,  as  more  elements  can  be  added  to  increase 
the  power  without  altering  the  operating  frequency.  It  is  also  free  of  the  multi-moding 
problem  which  plagues  traditional  combining  techniques,  because  each  element  is  designed 
to  operate  at  only  one  frequency.  Since  the  oscillators  are  electrically  separated  from  their 
neighbors,  they  can  be  designed  for  maximum  output  power  and  efficiency.  There  is  also 
more  flexibility  in  the  element  spacing  than  with  the  grid  approach,  so  that  higher  effective 
radiated  powers  per  element  are  possible.  Finally,  the  coupled  oscillator  architecture  exhibits 
graceful  degradation,  meaning  that  if  one  oscillator  fails  the  array  will  continue  to  operate 
normally,  albeit  with  slightly  deteriorated  performance.  This  is  in  contrast  to  the  oscillator 
grids,  where  failure  of  one  device  can  lead  to  failure  of  an  entire  row  or  more. 
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2.  Summary  of  the  Most  Important  Results 
2.1  Active  Radiating  Elements 

The  building  blocks  for  the  coupled  oscillator  arrays  are  the  individual  oscillator  elements, 
consisting  of  an  active  device  (Gunn,  IMPATT,  RTD,  FET,  etc.)  integrated  directly  into  a 
radiating  element.  Several  novel  architectures  have  appeared  in  the  literature  [12-16]  which 
creatively  incorporate  an  active  device  in  a  planar  microstrip  antenna.  Important  figures 
of  merit  are  the  output  power,  efficiency,  and  packing  density.  The  experimental  arrays 
described  later  were  constructed  using  two  different  active  microstrip  patch  antenna  designs, 
one  with  a  Gunn  diode  [12]  and  the  other  with  a  MESFET  [16],  as  shown  in  figure  3.  The 
patch  antenna  is  a  useful  structure  for  this  purpose,  since  the  devices  can  be  integrated 
vertically  and  heatsinking  is  facilitated  by  the  ground  plane. 


Figure  3  — Active  patch  antennas  using  two-  and  three-terminal  devices,  (a)  A  Gunn 
or  Impatt  diode  is  mounted  between  the  ground  plane  and  the  patch,  (b)  An  FET  is 
mounted  across  the  narrow  gap,  with  source  leads  grounded  through  the  substrate.  Bias 
lines  are  also  shown.  These  elements  are  simple  to  design  and  easy  to  fabricate,  making 
them  attractive  for  use  in  large  arrays. 


The  Gunn/patch  element  was  constructed  using  the  topology  of  figure  3a.  The  device  is 
located  at  the  point  where  it's  impedance  is  matched  to  that  of  the  patch.  This  position 
can  be  found  using  a  first-principles  time-domain  simulation  of  the  device  [IT]  and  a  suitable 
model  for  the  patch  input  impedance  (such  as  the  cavity  model).  Alternatively  a  semi- 
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Angle,  degrees 

Figure  4  —  E-  and  H-plane  radiation  patterns  for  a  7.6  GHz  active  patch  using  a  Gunn 
diode.  The  patterns  agree  with  expected  patch  behaviour,  although  the  H-plane  cross¬ 
polarization  is  abnormally  large.  This  could  be  due  to  the  large  size  of  the  packaged 
diode. 


empirical  approach  can  be  used  [18],  where  an  approximate  impedance  is  assumed  for  the 
device  and  the  proper  location  is  found  using  this  impedance.  The  semi-empirical  method  is 
convenient  for  hybrid  X-band  circuits,  but  the  more  involved  computer  modeling  is  important 
for  millimeter-devices  which  must  integrated  monolithically  into  the  antenna  and  hence  are 
more  difficult  and  expensive  to  make.  Bias  is  applied  at  a  low  impedance  point  on  the  patch. 
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A  FET/patch  element  was  also  developed  for  power-combining  arrays  [16],  and  is  shown 
in  figure  3b.  This  structure  is  not  a  conventional  patch  antenna,  since  it  uses  two  low- 
impedance  microstrip  lines  coupled  by  a  narrow  gap.  The  device  is  mounted  in  a  common- 
source  configuration  across  this  gap.  with  the  source  leads  grounded  through  the  substrate. 
The  feedback  capacitance  from  the  gap  is  sufficient  to  make  the  device  unstable  (and  hence 
useful  as  an  oscillator),  while  the  open  circuited  lines  provide  a  good  conjugate  match  to  the 
load  (radiation  resistance).  The  FET/patch  element  has  a  demonstrated  efficiency  of  26%, 
but  otherwise  behaves  similarly  to  the  Gunn/patch  element. 

Typical  radiation  patterns  for  a  Gunn/patch  element  are  shown  in  figure  4.  This  element 
measured  0.45  by  0.65  inches,  with  the  diode  0.15"  from  the  edge,  and  was  fabricated  on 
60  mil,  er  —  4.1  substrate.  Commercially  available  low-power  MA/COM  packaged  Gunn 
diodes  (MA  49104)  were  used  in  this  work.  Typical  bias  was  12  volts  @250  mA,  with 
an  overall  DC-to-RF  efficiency  of  1%.  Linear  polarization  and  measured  antenna  patterns 
were  consistent  with  expected  patch  behavior.  The  high  cross-polarization  in  the  H-plane 
measurement  indicates  excessive  coupling  to  other  patch  modes,  which  is  attributed  to  the 
large  size  of  the  diode  package  in  relation  to  the  antenna  —  this  effect  should  be  reduced  in 
a  monolithically  fabricated  array  [12],  since  no  bulky  package  would  be  used, 


Element  Separation,  x/X 

Figure  5  —  Oscillator  coupling  measurement  for  a  Gunn/patch  in  the  E-plane.  The 
theory  curve  is  calculated  as  described  in  the  text. 
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The  tuning  curves  (frequency  and  output  power  versus  bias  voltage)  for  the  circuit  were 
very  similar  to  other  published  measurements.  Bias  tuning  can  be  used  to  compensate  for 
small  discrepancies  between  the  elements,  which  is  important  because  the  proper  operation 
of  the  array  requires  that  the  elements  have  nearly  identical  characteristics.  For  an  array  of 
elements,  a  partial  reflector  (such  as  a  dielectric  slab)  positioned  above  the  array  facilitates 
the  mutual  injection-locking  of  the  devices  and  helps  establish  the  proper  phase  relationships 
between  the  elements.  Figure  5  shows  the  measurements  of  a  single  Gunn/patch  element 
when  a  dielectric  slab  is  moved  above  the  device. 

2.2  Coupled  Oscillator  Theory 

As  part  of  the  design  of  arrays  an  elegant  treatment  of  oscillator  coupling  was  developed 
based  on  Adler’s  equation  for  injection- locking  [19].  An  anlalysis  of  systems  using  using 
coupled  or  “inter-injection-locked”  microwave  oscillators  has  been  previously  published  [19, 
20].  These  analyses  are  elegant  but  are  impractical  for  power-combining  arrays  containing 
hundreds  or  thousands  of  elements.  The  theory  gave  us  considerable  physical  insight  into 
the  design  and  operation  of  these  arrays.  This  approach  was  motivated  by  recent  work  in 
low-frequency  coupled  oscillators  [24,25].  In  that  work,  an  unusual  mathematical  model  was 
postulated  for  the  oscillator,  which  dissociates  amplitude  and  phase  dynamics.  It  was  then 
argued  that  the  steady-state  behavior  is  a  function  of  the  phase  dynamics  alone,  and  after 
introducing  a  suitable  coupling  term,  a  compact  analytical  result  for  coupled  oscillators  was 
derived.  However,  the  models  were  chosen  without  any  correspondence  to  physical  reality, 
and  the  coupling  mechanism  was  assumed  instantaneous,  which  is  not  a  valid  assumption  at 
high  frequencies.  A  summary  of  the  theory  is  given  in  Appendix  I. 

The  theory  tell  us  that  stable  solutions  require  that  the  coupling  signals  between  elements 
are  sufficiently  strong  to  lie  in  the  locking  range  of  the  oscillators,  and  the  coupling  phase  is 
a  multiple  of  2tt.  This  theory  has  been  confirmed  by  a  series  of  array  measurements. 

2.3  Gunn  Diode  Arrays 

The  first  experimental  array  of  weakly  coupled  oscillators  was  a  16-element  array  using  Gunn 
diodes,  shown  in  figure  6  [10,  12].  The  individual  elements  were  designed  semi-empirically, 
as  described  previously.  This  array  design  uses  individual  bias  to  each  device,  which  was 
required  due  to  device  non-uniformities.  Elements  of  the  array  are  spaced  a  half  free  space 
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wavelength  apart,  a  distance  which  was  initially  selected  based  on  curves  in  (3j.  The  array 
was  designed  on  a  60  mil  substrate  with  er  =  4.1. 


Top  View 


Patch  with 
Diode 

Bias  Lines 


Brass 

Heatsink  and 
Groundplane 


Dielectric 

Substrate 


Figure  6  —  Diagram  of  Gunn  diodes  mounted  into  a  4x4  array  of  microstrip  patches. 
The  brass  block  serves  as  a  groundplane,  heatsink,  and  DC  bias  return.  Individual  bias 
to  each  element  is  applied  at  an  RF  null. 


Each  diode  was  first  biased,  one  at  a  time,  to  establish  a  common  operating  frequency.  These 
individual  biases  were  then  applied  simultaneously.  Single  frequency  operation  was  verified 
with  a  spectrum  analyzer,  as  shown  in  figure  7b.  Spectra  resembling  figure  7a  result  when 
the  elements  are  not  all  in  synchronization  (this  effect  has  been  fully  explained  using  Adler’s 
equation  [22,23]).  Slight  differences  in  diode  characteristics  and  diode  placements  made  the 
simultaneous  injection  locking  a  delicate  operation.  As  previously  suggested,  the  addition  of 
a  1  inch  thick  dielectric  slab  ( er  —  4.0)  above  the  array  facilitated  the  injection-locking. 

The  individual  elements  exhibited  a  frequency  tuning  from  9.5  to  10.0  GHz  versus  bias  volt¬ 
age.  The  best  results  for  the  array  were  obtained  at  a  frequency  of  9.6  GHz,  which  is  within 
4%  of  the  design  frequency.  Sharp  patterns  in  both  the  E-  and  H-plane  corresponding  to  a 
directivity  of  17  dB  were  measured  at  this  frequency  (figure  8),  and  indicate  in-phase  oper¬ 
ation.  Some  frequency  tuning  was  observed  by  adjusting  the  position  of  the  dielectric  slab, 
but  this  effect  is  limited  to  the  maximum  locking  range  of  the  array  («  100  MHz)  and  also 
changes  the  radiation  patterns  substantially. 

A  maximum  received  power  of  9.6  mW  was  obtained  at  a  distance  of  1.1  meters  away  from 
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Figure  7  —  (a)  Spectrum  of  the  Gunn  array  on  the  threshold  of  synchronization,  and 
(b)  at  full  synchronization.  These  measurements  were  made  with  the  dielectric  reflector 
in  place. 


the  array,  using  a  l9dB  pyramidal  horn.  The  total  radiated  power,  estimated  from  the 
measured  E-  and  H-pattems,  was  415  mW.  The  method  of  estimating  total  power  from 
the  principal  radiation  patterns  is  admittedly  inaccurate,  but  based  on  comparisons  with 
theoretical  models  we  believe  that  the  figure  presented  here  may  actually  underestimate  the 
total  power  by  as  much  as  16%.  This  total  power  figure  gives  26  mW  per  device,  which 
is  consistent  with  the  25  mW  rating  of  the  diodes  —  considerably  more  power  could  be 
obtained  using  higher-power  devices.  The  above  data  correspi  >nds  to  an  EIRP  of  22  Watts. 
The  overall  DC  to  RF  conversion  efficiency  was  low,  typical  of  Gunn  diodes,  around  1%. 


2.4  MESFET  Arrays 

While  t.ie  Gunn  diode  array  has  demonstrated  the  concept  of  power-combining  using  the 
weakly-coupled  oscillators,  the  low  efficiency  would  be  a  major  disadvantage  in  many  ap¬ 
plications.  Much  higher  efficiency,  larger  tuning  range,  and  better  noise  properties  can  be 
obtained  using  FET  devices.  Furthermore,  the  use  of  FETs  also  invites  optical  control 
control  of  the  array  [26]- 

The  experimental  MESFET  array  is  depicted  in  figure  9.  The  individual  element  design  was 
previously  discussed,  and  each  uses  a  genera)  purpose  Fujitsu  device  (FSX02).  Half  as  many 
bias  lines  were  used  in  this  array  design,  with  bias  isolation  between  the  elements  provided 
using  a  6-turn  coil.  The  gate  resistor  was  found  neccessary  to  supress  bias  circuit  oscillations. 
The  array  was  fabricated  on  a  93  mil  Duroid  5870  substrate  (cr  =  2.33). 
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Angle,  degrees 

Figure  8  —  Measured  E-  and  H-plane  patterns  at  9.6  GHz  for  the  active  array  of 
figure  6.  The  theoretical  results  are  calculated  by  combining  the  pattern  of  a  single 
patch  with  a  4x4  array  factor.  The  dielectric  slab  above  the  array  (er  =  4)  has  a  small 
effect  on  the  patterns.  Good  qualitative  agreement  between  the  measured  and  calculated 
curves  indicates  that  the  elements  are  nearly  in  phase  with  similar  amplitudes. 
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As  with  the  Gunn  array,  the  power-on  sequence  was  to  tune  each  group  of  elements  individu¬ 
ally  to  set  a  common  operating  frequency,  and  then  apply  DC  power  to  all  elements  at  once. 
Varying  the  reflector  element  spacing  is  then  usually  sufficient  to  enforce  mutual  synchroniza¬ 
tion.  With  single  frequency  operation  verified  at  8.27 GHz  (within  3%  of  the  8  GHz  design 
frequency  )  using  a  spectrum  analyzer,  the  patterns  were  then  measured.  These  patterns 
were  found  to  be  similar  to  the  Gunn  diode  array  and  closely  correspond  to  the  expected 
pattern  when  the  element  are  all  in  phase.  Linear  polarization  measurements  indicate  low 
cross-polarization  levels.  A  maximum  received  power  of  1.4  mW  at  a  distance  of  2.23m  was 
measured  using  a  19 dB  pyramidal  horn,  yielding  an  EIHP  of  10  W.  A  total  radiated  power 
of  184  mW  was  estimated  from  the  pattern  measurements,  giving  11.5  mW  per  device  at  26% 
efficiency.  The  corresponding  directivity  is  I7.2dB.  The  results  of  the  MESFET  array  are 
seen  to  be  similar  to  that  of  the  Gunn  array,  but  with  a  much  higher  DC-to-RF  efficiency, 
which  makes  it  a  more  attractive  design.  Better  efficiency  and  somewhat  higher  output 
power  per  device  could  be  obtained  by  using  gate  bias  to  each  element. 

The  MESFET  array  exhibited  a  broader  tuning  range  than  the  Gunn  array,  approximately 
400  MHz  with  varying  reflector  position.  Much  quieter  operation  was  also  observed,  which 
is  characteristic  of  FET’s.  The  spectrum  of  the  MESFET  array  is  shown  in  figure  10,  along 
with  the  spectrum  when  the  array  is  externally-locked  to  an  HP835QB  sweeper  (83592B 


-0.5  -0.4  -0.3  -0.2  -0.1  0  0.1  0.2  0.3  0.4  0.5 

Frequency  Offset.  MHz 

Figure  10  —  FYee-running  and  externally- locked  spectra  of  the  MESFET  array.  The 
locking  signal  was  injected  through  a  pyramidal  horn,  which  illuminated  the  array  from 
a  distance  of  1  meter. 

plug-in).  Phase  noise  of  -78 dB  /Hz  at  100 kHz  for  the  free-running  case  was  measured  using 
a  spectrum  analyzer. 

2.5  Mode-locked  Arrays 

A  mode-locked  system  is  created  when  a  number  of  different  equally-spaced  frequencies, 
or  modes,  are  simultaneously  produced  and  locked  in  phase  with  each  other  [27).  Duing 
the  course  of  this  contract,  we  discoverd  that  the  weakly  coupled  arrays  could  be  operated 
in  an  interesting  mode-locked  configuration.  This  mode-locking  can  be  accomplished  by 
arranging  the  frequencies  of  each  oscillator  in  an  array  to  produce  a  spectrum  with  equally 
spaced  components,  as  shown  in  figure  11.  The  mutual  coupling  between  the  array  elements 
then  serves  to  lock  the  frequencies  and  phase-relationships  together.  This  is  a  ‘passive' 
mode-locking  scheme,  but  ‘active’  mode-locking  may  be  possible  for  increased  stability  and 
reliable  operation. 

Recent  experimentation  has  revealed  an  additional  feature  of  this  system,  namely  that  the 
pulse  train  is  effectively  scanned  above  the  array  as  in  a  phased  array  source.  Although 
previously  unexpected,  this  phenomenon  follows  directly  from  the  mode-locking  theory  and 
elementary  array  theory.  The  scan  coverage  is  determined  by  the  element  spacing,  while 


16 


H  Am  !■ — 


Figure  11  —  (a)  Frequency  spectrum  of  a  5-eiement  mode-locked  oscillator  array,  (b) 
Phasor  description  in  the  complex  plane  at  one  instant  of  time. 


the  pulse  repetition  frequency  is  determined  by  the  spacing  of  the  Spectral  components.  An 
oscillator  array  operated  in  this  mode  could  potentially  function  as  the  transmitter  in  a 
phased-array  radar  application,  without  the  need  for  any  phase  shifters.  In  addition,  since 
the  individual  elements  are  operating  in  a  CVV  mode,  there  are  no  thermal  restrictions  on 
the  pulse  repetition  frequencies  and  no  need  for  injection-priming  of  the  array. 


Figure  12  —  Experimental  setup  for  measuring  the  scanned  pulse.  A  fixed  receiving 
horn  at  the  broadside  position  collects  the  output  signal  for  observation  on  a  spectrum 
analyzer  and  oscilloscope  while  the  array  is  rotated. 


The  total  electric  field  from  an  array  of  M  =  2m  +  1  oscillators  spaced  at  frequencies 
u  ~  uq  -  iAu  with  locked  phases  such  that  <&  -  <k_i  =  8<j>,  (t  =  -m, . . . ,  m)  can  be  written 


17 


as  [27] 

E(t)  =  £g 

where  £'  is  the  elapsed  time  from  the  point  at  which  all  the  oscillators  have  the  same  phase, 
and  Eq  is  the  amplitude,  assumed  identical,  of  each  of  the  oscillators.  Equation  (1)  represents 
a  carrier  at  frequency  a?o  modulated  by  a  train  of  pulses.  The  peak  power  of  each  pulse  is 
EqM2,  so  that  an  array  of  10  elements  would  produce  a  peak  power  of  100  times  that  of  a 
single  element. 

Preliminary  measurements  using  a  5-element  linear  array  at  10 GHz  show  promising  results 
[28,29].  The  experimental  setup  is  illustrated  in  figure  12.  A  receiving  horn  is  located  at  the 
broadside  position,  while  a  time-reference  trigger  signal  at  the  pulse  repetition  frequency  is 
derived  by  capacitively  coupling  to  the  bias  lines  on  the  array.  The  output  signal  from  the 
horn  is  observed  on  a  spectrum  analyzer  and  a  high-speed  oscilloscope  as  the  array  is  rotated 
in  the  scanning  plane. 
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Figure  13  —  Spectrum  of  three  mode-locked  microwave  oscillators.  Additional  side¬ 
bands  arise  from  frequency  pulling  effects  [28]. 

An  experimental  observation  of  both  the  mode-locked  pulse  phenomena  and  the  scanning 
was  made  by  adjusting  the  oscillation  frequencies  of  each  element  (using  the  bias  supplies) 
to  give  equally  spaced  spectral  components.  This  step  was  facilitated  using  the  spectrum 
analyzer  as  shown  in  figure  13.  Note  from  figure  13  that  additional  sidebands  are  produced 
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as  a  result  of  the  mutual  pulling  effects.  The  pulse  repetition  rate  is  inversely  proportional 
to  the  frequency  difference,  which  was  constrained  by  the  locking  range  of  the  oscillators  to 
50  MHz  or  more. 
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Figure  14  —  Measured  (•)  and  theoretical  (solid  line)  time  dependence  of  the  output 
signal  power  envelope  of  five  mode-locked  oscillators. 

Once  the  proper  array  spectrum  was  verified  on  the  spectrum  analyzer,  a  mode-locked  pulse 
was  dearly  observed  on  the  oscilloscope  as  shown  in  figure  14.  In  this  figure,  only  the 
envelope  of  the  output  signal  is  shown  for  clarity.  The  theory  curve  presented  in  figure  14  is 
calculated  from  equation  (1),  and  shows  excellent  agreement  to  the  measurements. 

When  the  array  is  rotated  (as  in  figure  12)  the  pulse  is  observed  to  shift  in  time,  corresponding 
to  the  pulse  being  scanned  above  the  array.  This  is  shown  in  figure  16  for  three  different 
scan  angles.  At  large  angles  the  pulse  strength  diminishes  due  to  the  radiation  pattern  of 
the  individual  patch  antennas  in  the  array. 
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Figure  15  —  Pulse  waveforms  for  a  three  element  active  patch  array  at  three  scan 
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Appendices 


Appendix  1  Coupled  Oscillator  Theory 

The  assuption  that  individual  oscillators  in  an  array  are  only  slightly  perturbed  from  their 
free- running  state  by  the  presence  of  the  other  oscillators  leads  naturally  to  Adler’s  equation 
for  injection  locking  [21-23], 


doo  Ain1  •  /  i  ,  \  /  , 

-fa  =  — ^2g  sin(^  ~  +  ~  ^ 


(A.l) 


where  <fio  —  phase  of  oscillator.  cinj  =  phase  of  injected  signal,  uq  =  free-running  frequency  of 
oscillator,  W;nj  =  frequency  of  inje  -ed  signal  from  an  external  or  neighboring  oscillator,  Ao  = 
free-running  amplitude  of  oscillator  (voltage),  Ainj  =  amplitude  of  injected  signal  (voltage), 
and  Q  =  the  external  Q  of  the  oscillator  circuit.  The  phase  of  the  oscillator  is  defined  relative 
to  the  injected  signal  frequency,  so  that  the  steady-state  solution,  d<t>o/dt  —  0,  corresponds 
to  the  oscillator  being  frequency- locked  to  the  injected  signal.  Adding  the  injected  frequency 
to  both  sides  of  (A.l)  and  noting  that  uj  —  4-  dfo/dt  —  the  instantaneous  frequency  of 

the  oscillator,  yields 

w  =  w o  1  -  sin(0o  -  (A2) 

For  mutually  syncronized  arrays,  (A. 2)  must  be  modified  to  include  the  effects  of  coupling 
delay  and  several  simultaneous  injected  signals.  First,  the  coupling  between  individual  oscil¬ 
lators  can  be  expressed  in  terms  of  a  complex  coupling  coefficient.  This  is  just  the  scattering 
parameter  Sij  of  the  circuit  connecting  the  two  oscillators  i  and  j  (reciprocity  is  assumed, 
Sjj  =  Sji).  The  magnitude  and  phase  of  this  coupling  coefficient  are  written  separately  as 
=  Aij  exp(jd>tJ).  Secondly,  superposition  is  used  to  account  for  several  simultaneously 
injected  signals.  Incorporating  the  above  for  a  system  of  N  oscillators  gives  a  modified  form 
of  Adler’s  equation  for  the  ith  oscillator, 


w  =  1  sin(0i-0j  +*ii) 

L  ZQ'  /i‘ 


i  =  1,2 . JV 


(A.3) 


where  the  instantaneous  frequency  u  of  oscillator  i  is  the  same  for  all  oscillators  in  the 
system  at  synchronization.  Equation  (A. 3)  describes  a  set  of  equations  which  can  be  used 
to  determine  the  steady-state  of  the  system.  If  the  free-running  parameters  are  given,  the 
unknowns  are  the  phase  differences  and  operating  frequency  of  the  array.  Conversely  one 
can  solve  for  the  frequencies  and  amplitudes  which  give  a  prescribed  phase  distribution. 
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Some  simple  analytical  results  can  be  derived  from  (A. 3)  by  considering  a  linear  chain  of 
oscillators  with  only  nearest  neighbors  interactions  (the  more  general  situation  of  a  two- 
dimensional  array  is  similar  but  notationally  more  complicated).  Assuming  that  the  coupling 
is  the  same  between  adjacent  elements  in  the  array,  Ay  =  A  and  ^l}  =  4>.  Furthermore,  let 
Q  —  Qi,  A'  =  A /2Q,  pi  =  and  Arii  =  —  4>i~ i-  The  set  of  governing  equations 

becomes 


tjj  —  [1  —  A'piSin(<f>  -r  AO,)  —  A'3in(3>  —  Ad>i-*-i)/P»+i]  n=l,2,...,jV  (A. 4) 


Note  that  p\  —  1  /pn+\  —  0-  This  is  the  same  form  obtained  in  [25],  except  for  the  presence  of 
the  coupling  phase  term.  However  this  additional  term  renders  the  equations  unsolvable  by 
conventional  linear  techniques.  To  solve  these  equations,  further  assumptions  must  be  made 
such  as  a  small  angle  restriction  where  sin(Ad>i)  «  A <f>it  or  a  nonlinear  root-finding  procedure 
could  be  used.  For  small  numbers  of  elements,  the  equations  can  be  solved  analytically. 


In  general  there  are  many  possible  phase  distributions  which  satisfy  this  set  of  equations,  but 
not  all  are  stable  solutions.  Stability  of  these  modes  wm  investigated  using  a  perturbation 
analysis.  Recalling  that  cj  =  4-  dfii/dt  a  new  differential  equation  is  obtained, 

^A 4>i  =  Ui  [1  -  A V,  sin(<f>  4-  A &)  -  A'sin(4>  -  Ao,+i)/pi+il  ^  ^ 

-  Ui- 1  [1  -  A'pi_isin(^  4-  A0i_i)  -  A'sin(4>  -  A d>i)jpi) 


The  phase  is  then  perturbed  by  a  small  amount  by  letting  Ao, 
algebra  this  leads  to 


A0»  4-  After  some 


—6i  =  Uj6i_i  4-  bi6i  4-  Ci6i+i 


i  =  2, 3, . . . ,  JV 
where  a*  =  A'pi-M--!  cos(^>  4-  Ad>»-i) 

bi  =  —\,UipiCos(’&  4-  A 0i)  —  cos(<$  —  Ad h)/p% 


(A.6) 


Cj  =  A'u/iCos(^  -  A <t>i+i)/pi+l 

This  is  of  the  form  of  a  matrix  equation,  dSjdt  =  AS,  where  A  is  a  tridiagonal  matrix,  and 
6  is  a  column  vector  with  elements  Si.  A  stable  solution  for  the  phase  distribution  requires  S 
to  decay  with  time.  This  is  satisfied  when  all  the  eigenvalues  of  A  have  negative  real  parts. 
In  all  of  the  numerical  solutions  of  (A. 4)  and  (A. 6)  which  have  been  carried  out  so  far,  the 
stability  criteria  (A. 6)  has  eliminated  all  but  one  or  two  of  the  possible  modes  given  by  (A. 4). 
This  observation  has  resisted  an  analytical  justification  so  far.  It  has  also  been  observed  that 
for  certain  values  of  the  coupling  phase  there  are  no  stable  modes  of  operation. 
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To  illustrate  the  theory  and  derive  some  important  qualitative  results,  consider  a  simple 
chain  of  four  similar  oscillators.  From  (A.4), 

u  =  u\  [1  -  A'  sin(#  -  A&)] 
u  =  v<2  [1  -  A'sin(#  +  A#j)  -  A'sin(#  —  A03)] 

(/IT) 

u  =  a/3  [1  —  A'  sin(#  +  A03)  —  A'  sin(#  -  A^)] 
u  =  u?3  [i  —  A'  sin(#  +  A^4)j 

assuming  that  the  oscillator  amplitudes  are  the  approximately  the  same  so  that  pt  ~  1.  This 
is  usually  a  good  approximation  in  practice.  From  (A.  7),  if  the  amplitudes  are  similar,  then 
the  final  oscillation  frequency  of  the  array  will  lie  within  the  locking  range  of  each  oscillator; 
that  is,  the  locking  ranges  of  every  oscillator  must  overlap  to  some  extent.  Since  these 
are  weakly  coupled  oscillators,  the  locking  bandwidth  is  small,  hence  the  elements  must  be 
nearly  identical  to  satisfy  (A. 7).  It  is  important  to  note,  from  (A.2),  that  since  the  oscillator 
phase  undergoes  a  180°  variation  over  the  locking  range,  small  differences  in  the  free-running 
frequencies  can  have  a  large  effect  on  the  solution.  Assuming  perfectly  identical  oscillators 
{ ui  =  u-2  =  u3  =  un),  gives 

sin(#  —  A^)  =  sin(#  +  A<&)  +  sin(#  —  A<A?) 
sin(#  +  A<fo)  4-  sin(#  -  A<£3)  =  sin(#  +  A<fo)  +  sin(#  —  A<£4)  (A8) 

sin(#  +  A <f>3)  +  sin(#  —  A<^4)  =  sin(#  -f  A<Aa) 

This  shows  that  end  effects  will  be  important,  at  least  for  small  arrays,  since  the  elements 

on  the  periphery  see  quite  different  injection  signals.  There  are  many  possible  solutions 

to  (A.8);  for  phase-coherent  power-combining,  of  most  interest  is  the  solution  for  which 
Ad>2  =  A^>3  =  Ad>4  =  0.  For  this  solution,  (8)  is  satisfied  when  #  =  nn,  where  n  =  0, 1, 2  — 
From  the  stability  analysis,  the  following  eigenvalue  equation  is  obtained 

— 2  cos#  —  e  cos#  0 

cos  #  — 2cos#  — e  cos  #  —0  (A. 9) 

0  cos#  —  2  cos#  —  e 

The  eigenvalues  e  can  have  negative  real  parts  only  if  cos#  >  0,  so  a  stable,  in-phase 
mode  is  only  possible  if  #  =  0, 27r,....  If  free-sp'-’e  coupling  is  predominant,  this  means 
that  the  elements  must  be  spaced  at  mutiples  of  one  wavelength.  Such  spacing  is  generally 
unacceptable  because  of  grating  lobes  in  the  antenna  patterns.  Furthermore,  it  has  been 
shown  that  close  spacing  is  desirable  for  large  arrays  for  efficient  power  combining  [3J. 
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Appendix  2.  Mode-Locked  Array  Theory 

For  an  array  of  resonant  radiating  structures,  a  very  simple  circuit  model  has  been  found  to 
adequately  predict  experimentally  observed  effects.  The  model  consists  of  an  active  device, 
such  as  a  Gunn  diode,  embedded  in  a  series  resonant  circuit.  It  is  assumed  that  the  active 
device  can  be  modelled  by  an  effective  negative  resistance,  whose  magnitude  is  amplitude- 
dependent  but  frequency-independent  to  first  order.  A  schematic  for  this  circuit  is  shown  in 
figure  5,  where  an  additional  source  is  included  for  externally  injected  signals. 
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Figure  16  —  Simple  circuit  model  for  a  single  oscillator.  The  device  is  represented  by 
a  lumped  negative  resistance.  A  single-tuned  resonant  circuit  is  assumed,  so  that  har¬ 
monic  effects  will  be  minimal.  The  source  Vinj  accounts  for  interactions  with  neighboring 
oscillators. 


The  circuit  equation  corresponding  to  figure  16  is 


,nj  Rl  dt  + 


V0  dt  +  V0 


(A.10) 


where  V0  is  the  output  voltage  amplitude  and  Vinj  is  the  injected  signal  due  to  neighboring 
oscillators.  For  convenience,  let  <5(Vq)  =  1  —  Rj/Rl,  where  the  dependence  of  the  negative 
resistance  on  amplitude  has  been  explicitly  written.  If  Vinj  =  0,  then  the  oscillation  condition 
would  require  <5(Vo)  =  0,  however  when  there  is  an  injected  signal  present  this  parameter  is 
non-zero  to  accommodate  the  additional  energy  being  supplied  to  the  circuit.  For  a  series 
ti  ned  circuit 


2 1  Rl  _  tuo 

U°~~LC  T  ”  ~Q 


where  Q  is  the  quality  factor  of  the  passive  resonating  circuit.  The  Q-factor  is  assumed  to 
be  large  enough  ( Q  >  10)  so  that  the  output  signal  will  be  close  to  the  natural  resonant 
frequency  of  the  circuit, 


V0  =  A0(t)  exp  {j  (uot  +  <M0)}  =  A0cJu^t  (A. 11) 
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Following  van  der  Pol’s  approach  (or  the  method  of  averages),  the  amplitude  and  phase 
variables  are  allowed  to  be  slowly  varying  functions  of  time,  leading  to  the  set  of  equations 

oRe{^}  (A,2„) 

f  =  (A12W 

Equation  (A.  12b)  is  recognized  as  Adler’s  equation  for  injection-locking. 

The  negative  resistance  appearing  in  (A.  12),  via  the  term  <5(Ao),  is  a  time-averaged  value 
described  by 


where  V'd  and  /d  are  the  voltage  and  current  at  the  terminals  of  the  device.  Consistent 
with  the  implicit  assumptions  of  our  circuit  model,  the  simplest  nonlinearity  which  yields 
a  sinusoidal  oscillation  is  a  cubic  relationship  between  the  current  and  voltage,  which  gives 
an  approximate  expression  for  the  amplitude  dependence  of  negative  resistance  as  Rj{A)  = 
a  —  bA 2.  The  amplitude  dependence  of  negative  resistance  is  the  mechanism  responsible 
for  the  amplitude-limiting  of  any  real  oscillator.  Normalizing  this  expression  such  that  the 
free-running  oscillations  occur  when  do  =  1  gives 

6(/t0)  =  -m(1  -  V)  (A  13) 

where  n  =  6/ftz,  is  a  parameter  representing  the  oscillator.  With  this  relation,  (A.  12)  is  seen 
to  be  closely  related  to  the  van  der  Pol  equation.  In  the  following  discussions  (A.  13)  will  be 
used  to  represent  the  amplitude  dependance,  and  n  is  taken  as  an  empirically  determined 
quantity. 

For  a  coupled  system  of  such  oscillators,  the  coupling  between  elements  can  be  described  by 
a  complex  coupling  coefficient,  where 


Kij  =  Aijexp(-jrl>jj) 


for  the  coupling  between  oscillators  i  and  j  in  the  system.  Reciprocity  is  assumed  in  the 
following  discussion  =  Kji).  The  injected  signal  seen  by  the  ith  oscillator  is  then 

/v 


;=i 

i  t4* 
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where  N  is  the  number  of  oscillators  in  the  system  and  Vj  is  the  output  signal  of  the  yth 
oscillator.  With  A  and  u>i  denoting  the  amplitude  and  free-running  frequency  of  the  ith 
oscillator,  and  using  the  instantaneous  phase  =  ujit  -f  fa,  a  system  of  equations  describing 
the  coupled  oscillator  dynamics  can  be  derived,  giving 

,v 

M  =  -f ^  «•(♦«  +  *•- »j) 

i  =  l, 2,. ..,1V  (A14) 

+  &  -  0j) 

Experiments  have  indicated  that  nearest-neignbor  interactions  among  oscillators  in  a  ra- 
diatively  coupled  system  are  the  most  important,  as  depicted  in  figure  17.  Assuming  the 
coupling  to  be  the  same  between  adjacent  oscillators  in  the  system  yields  ^  and 

<Pij  =  <P.  Using  this  simplification,  and  writing  out  5(A)  explicitly,  simplifies  (A. 14)  to 


A  =  “Ml  ~  A2)A  +  ^  AjCos(<&  +  di-dj) 

i  i*' 

i+l  , 

j¥> 


i  =  1,2,.. 


(A.lo) 


The  relative  phase  differences  between  the  oscillators  is  of  most  interest,  so  let  Ad\  =  —  #,-i 

and  Ao>i  =  Wj—  Wj_i-  Furthermore,  note  that  the  equations  (A. 14)  are  only  valid  to  first  order 
in  l/Q,  and  for  Aw,  <§:  04,  Under  these  assumptions  the  quantity  <jJi/2Q  is  approximately 
the  same  for  every  oscillator.  Defining  9!  =  ujq/2Q  and  X1  =  u/0/2t?,  where  uo  is  an  average 
frequency  for  the  oscillator  ensemble,  yields 


A  =//(!  -  A3)A  +  A7<(A,A0)  i  =  \.2,...,N  (A  15a) 

Adi  =  Au \  4-  \'9i( A,  A 6)  i  =  2, 3, ... ,  .V  (A.156) 

where  /i(A,  Ad)  —  A-i  cos($  4-  Ad,)  +  A+i  cos($  —  A0i+1) 

g,(A,  AO)  =  sin($  4-  A#i_i)  4-  -  sin(<J>  -  A0,) 

A-i  A-i 

-  sin(<£  4-  Adi)  ~  ~  sin($  -  A0i+l) 

Ai  A* 

and  where  Aq  =  Av+i  =  0  is  assumed.  For  X  =  0  the  oscillators  are  uncoupled,  and  (A. 15) 
reduces  to  a  set  of  isolated  limit  cycle  oscillators  with  amplitudes  Ai  =  1  and  frequencies  a*. 
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Figure  17  —  A  chain  of  radiatively-coupled  oscillators,  with  bilateral  coupling  between 
nearest  neighbors. 


Equation  (A.  15b)  can  be  considered  a  modified  version  of  Adler’s  equation.  For  loose  cou¬ 
pling  and  a  narrow  distribution  of  oscillator  frequencies,  the  amplitudes  will  remain  close  to 
the  free-running  values,  and  the  osciilator  dynamics  will  be  essentially  governed  by  (A. lob). 
In  this  case  all  of  the  oscillators  in  the  system  may  lock  to  a  single  frequency,  associated 
with  stable,  fixed  points  of  (A. lob).  The  parameter  A'  in  (A. 15)  is  closely  related  to  the 
locking  bandwidth  of  the  oscillators.  Thus  for  widely  varying  oscillator  frequencies  such  that 
Aw,  >  ,V,  the  oscillators  may  not  lock  to  a  single  frequency.  In  this  case  the  amplitude 
dynamics  described  by  (A. 15a)  can  not  be  ignored,  and  both  A,  and  A0»  will  be  functions 
of  time  in  the  steady-state.  This  is  the  regime  in  which  mode-locking  can  occur. 
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